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ABSTRACT 

The functions are determined in a way that leads to the calculation of the a 
ctions of perfect recessing of reinforced concrete beams, with straight haunches, 
any kind of load and new features are presented for the CEINCI- 
LAB computer system 

with which is the moment and forces of perfect recessing in beams with straight 
haunches. 

Finally, two applications are made, one for static analysis of a flat gantry and one fo 
r spatial analysis of a structure with straight haunches. 

The two cases illustrate the use of CEINCI-LAB functions 

Keywords: Reinforced Concrete Straight Haunches. CEINCI-LAB computer system 


ACCIONES DE EMPOTRAMIENTO PERFECTO EN 
VIGAS CON CARTELAS LINEALES Y NUEVAS 
FUNCIONES DE CEINCI-LAB 

RESUMEN 

Se determinan las funciones de forma que conducen al calculo de las 
acciones de empotramiento perfecto de vigas de hormigon armado, con cartelas 
lineales, ante cualquier tipo de carga y se presentan nuevas funciones para el 
sistema de computacion CEINCI-LAB con las cuales se halla el momento y fuerzas 
de empotramiento perfecto en vigas con cartelas. 

Finalmente se realizan dos aplicaciones una para el analisis estatico de un 
portico plano y otra para el analisis espacial de estructuras con cartelas. En los dos 
casos se ilustra el uso de las funciones de CEINCI-LAB 

Palabras Claves: Cartelas Lineales de Hormigon Armado. Sistema de computacion 

CEINCI-LAB 
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1. INTRODUCCION 

El Edificio Administrativo de la Universidad de Fuerzas Armadas ESPE, es 
irregular en planta y elevacion, con grandes luces, por lo que fue disenado con vigas 
con cartelas lineales como la que se aprecia en la figura 1 a. Se aprecia que tiene 
vigas principales con cartela (que van de columna a columna) y vigas secundarias 
de seccion constante (en el centro de luz). Caiza (2005) 



(b) 


Figura 1 Estructuras de hormigon armado con cartelas; a) Edificio Administrativo de 
Universidad de Fuerzas Armadas ESPE; b) Bodega de Firmesa S.A., en 
Carcelen. 
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La bodega de la empresa Firmesa S.A., ubicada en Carcelen en Quito, tiene 
que cubrir grandes luces de 10 m., en los dos sentidos ortogonales. El dueno del 
Proyecto solo permitia tener columnas de hormigon armado en los porticos 
perimetrales. La figura 1 b, muestra un portico exterior perimetral con luces de 3 m, 
y 4 m, falta otra luz de 3 m, que no aparece en la fotografia; tambien se observan 
dos porticos interiores con vigas con cartelas lineales. Aguiar (2008). 
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Figura 2 Vista en planta de la bodega de la empresa Firmesa S.A., ubicada en 
Carcelen; b) Geometria de la viga con cartelas. 

De tal manera que los porticos exteriores si tienen columnas de hormigon 
armado pero los porticos interiores no tienen columnas, eso permite que puedan 
ingresar vehiculos a la planta baja con facilidad para dejar o llevar productos; en el 
segundo piso la idea inicial fue tambien tener una sola gran area. 

Para cubrir la luz de 10 m, en cada direccion se decidio tener vigas 
rectangulares de 40/40 cm, en la parte central de 4 m, y vigas con cartelas en los 
extremos de 3 m, de longitud; la dimension de las vigas en los extremos es de 40/60 
cm, ver figura 2. 

Para el analisis ante cargas verticales, de las vigas de hormigon armado 
con cartelas, se modelo como una malla espacial, la misma que se indica en la figura 
3, en la interseccion de las vigas no existe columnas. 
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Figura 3 Modelo numerico para analizar las vigas con cartelas de un Piso. 

En el Proyecto Estructural se calculo con las luces indicadas en la figura 2 
pero en la construccion se incremento ligeramente las dimensiones con unos 
voladizos. Esto es normal, de que en la construccion se modifique ligeramente el 
proyecto de ahf la necesidad del control y asesoramiento que debe dar el Proyectista 
Estructural en forma periodica. 

Lo cierto es que una buena alternativa estructural para cubrir grandes luces 
es que las vigas tengan cartelas. Ahora se debe decidir con que modelo de analisis 
se desea trabajar y por el momento se tiene dos casos desarrollados: i) Considerar 
a la viga con cartelas sin nudo intermedio como el que se muestra en la figura 4 a; 
ii) Considerar un nudo intermedio en la viga, de tal manera que se tienen dos 
elementos, el uno con cartela a la izquierda y el otro con cartela a la derecha, lo que 
se indica en la figura 4 b. 

En los dos modelos b el ancho constante de la viga; h la altura en el centro 
de luz; H la altura en el extremo de la viga; L v la longitud de la cartela; L la longitud 
total de la viga, para el modelo de la figura 4 a; L es la mitad de la longitud de la viga 
para el modelo de la figura 4 b. 

En Aguiar y Bonilla (2020), se presento la deduccion de las ecuaciones que 
conducen al calculo de la matriz de rigidez del elemento para los tres casos 
indicados en la figura 4, que son: i) Viga con dos cartelas; ii) Viga con cartela a la 
izquierda; iii) Viga con cartela a la derecha. 

Se presentaron tablas para quien quiera calcular manualmente los 
elementos de la matriz de rigidez en sistema de coordenadas 1 (sin considerar 
desplazamientos como cuerpo rfgido) y despues encontrar la matriz de rigidez en 
coordenadas locales y globales (Aguiar 2020). Se indico tambien las nuevas 
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funciones del sistema de computacion CEINCI-LAB que facilitan el calculo de la 
matriz de rigidez de una estructura. Bonilla (2019). 


b h 





Tipo = 0 Tipo = 1 

(b) 

Figura 4 Modelos de analisis de una viga de hormigon armado con cartelas; a) Sin 
considerar nudo en la mitad de la viga; b) Considerando nudo en la mitad 
de la viga, con lo que se tienen 2 elementos. 


Ahora, en este articulo, primero se presenta el marco teorico que conduce 
al calculo de las funciones de forma para encontrar las acciones de empotramiento 
perfecto para cualquier tipo de carga, para los dos modelos indicados en la figura 4 
b, que son: i) Viga con cartela a la izquierda; y, ii) Viga con cartela a la derecha. 


Despues se presentan tablas para hallar los coeficientes que llevan al 
calculo de las acciones de empotramiento perfecto, ante carga uniforme distribuida 
para el modelo de la figura 4 b. 


Finalmente se indican las nuevas funciones de CEINCI-LAB con las cuales 
se puede realizar el analisis estatico de un portico de hormigon armado con cartelas. 
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2. FUNCIONES DE FORMA 


De manera que se presente compatibilidad con CEINCI-LAB , el presente 
trabajo utiliza el modelo de la Figura 4 (b). Se adopta la nomenclatura mencionada 
en la misma y se considera los siguientes datos: 


h 


H 


P 


Lv 

~L 


Se puede encontrar la deduccion de ecuaciones de inercia en Aguiar (2020), 
Bonilla (2019), de modo que se tiene. 




Para: 0 < x < fiL 

/(x) — Ia* 



1 — a 
p * L 


x ) 


3 


• Para: flL < x < L 


/(x) — a 3 * Ia 

Partiendo de una carga Po=0 se deduce una ecuacion general de 
momentos integrando la carga. 


dM 2 

d 2 x 


= Po = 0 


dM 

dx 


A 


M(x) = A * x + B 

Utilizando la ecuacion de resistencia de materiales se tiene: 

dv 2 M(x) 
d 2 x E * I(x) 


• Para: 0 < x < pL 


d 2 v 

dx 2 


A * x + B 


E * Ia* 


1 — a 



( 1 ) 


Acciones de empotramiento perfecto en vigas con cartelas... 


241 


• Para: (2L < x < L 


d 2 v A * x + B 
dx 2 E * a 3 * Ia 

Para hallar las funciones de forma, es necesario resolver las ecuaciones 1 
y 2 por integrales sucesivas, de lo que se tiene: 

• Para: 0 < x < (3L 


6 > 0 ) 
= C - 


AEp* + ( BL 3 a - BL 3 )p 3 + x(2AL 3 a - 2AL 3 )(L 3 


(3) 


Ela((2L 2 a 2 - 4L 2 a + 2L 2 )(S 2 + x(4La 3 - 12 La 2 + 12 La - 4L)(S + x 2 (2a* - 8a 3 + 12a 2 - 8a + 2)) 


AL 3 (5 3 ln(Lp + x(a — 1)) 


AL 4 (3 4 + (BL 3 - BL 3 a)|3 3 


v(x) = ~- 


a 3 — 3a 2 + 3a — 1 (2La 3 — 6La 2 + 6La — 2L)P + x(2a 4 — 8a 3 + 12a 2 — 8a+ 2) 


E * Ia 


(4) 


+ Cx + D 


• Para: (3L < x < L 


Ax 2 - 


0(x) = 


+ Bx 


Elaa 3 


■+ E1 


(5) 


Ax 3 Bx 2 

v(x) = + E1 x + F 

Elaa 6 


( 6 ) 


Partiendo de estas ecuaciones, se debe utilizar varias condiciones de 
frontera, dependiendo de la reaccion que se desee encontrar. 

De manera que para encontrar el cortante en el nudo inicial, es necesario 
realizar un desplazamiento vl mientras el resto de las deformaciones son 0. 


x =0 0=0 

x = 0 v = vl 
x = L 9 = 0 
x = L v = 0 


Reemplazando las condiciones se obtienen las 4 primeras ecuaciones para 
el calculo de las constantes de integracion. 
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AL 4 /? 4 + (BL 3 a-BL 3 )/] 3 
E\a(2L 2 a 2 - 4L 2 a + 2L 2 )/? 2 


AL 3 f> 3 ln (Lf) AL*(1 4 + ( BL 3 - BL 3 a')fL 3 

a 3 — 3 a 2 + 3a — 1 + (2 La 3 — 6La 2 + 6 La — 2L)/? 
£3a 


( 7 ) 

( 8 ) 


AL 2 


+ £L 


£Iaa 3 


■ + E1 = 0 


( 9 ) 


AL 3 , BL 2 
6 + 2 
£Iaa 3 


+ EIL + F = 0 


( 10 ) 


Con el fin de completar el sistema de ecuaciones, se utilizan las condiciones 
que enlazan ambos tramos de la viga, estas son validas para encontrar todas las 
acciones de empotramiento perfecto. 


x — (3 * L 61 — 62 
x — (B * L vl = v2 


Reemplazando estas se tiene: 


L(2AL 3 a - 2AL 3 )/? 4 + +L 4 /? 4 + ( BL 3 a - BL 3 )p 3 

LIa(L 2 (2ft: 4 - 8c+ + 12a 2 - 8a + 2)/? 2 + L(ALa 3 - 12La 2 + 12La - 4L)/? z + (2L 2 a 2 - 4 L 2 a + 2L 2 )/? 2 ) 


+L 2 /? 2 

2 


+ BL/? 


Llaa 3 


+ E1 


( 11 ) 


AL 3 fl 3 ln(L(a - 1)/? + L/?) 


4L 4 /? 4 + (BL 3 - BL 3 a)p 3 


a 3 — 3a 2 + 3a — 1 L(2a 4 - 8a 3 + 12a 2 - 8a + 2)yg + (2La 3 - 6La 2 + 6La - 2L)/? 

Lla 


AL 3 (L 3 ^ BL 2 f> 2 

+ CL/? + D =-^^-+ EiL/? + F 


(12) 


Resolviendo el sistema de ecuaciones, se tiene las constantes de 
integracion para cortante en el nudo inicial. 
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A=(6*E*la*v1*(a-1 ) A 3*a A 3*(a A 2*p+a*p- 

2*p+2))/(L A 3*(6*a A 5*ln(a)*p A 4+6*a A 4*ln(a)*p A 4-12*a A 3*ln(a)*p A 4- 
14*a A 5*p A 4+22*a A 4*p A 4+a A 3*p A 4-13*a A 2*p A 4+5*a*p A 4-p A 4+12*a A 3*ln(a)*p A 3- 
22*a A 3*P A 3+36*a A 2*P A 3-18*a*P A 3+4*P A 3-6*a A 4*P A 2+24*a A 3*P A 2- 
36*a A 2*p A 2+24*a*p A 2-6*p A 2+2*a A 5*p-4*a A 4*p-4*a A 3*p+16*a A 2*p- 
14*a*p+4*p+a A 3-3*a A 2+3*a-1)) 

B=-(6*E*la*v1*(a-1 ) A 3*a A 3*(a*p A 2-p A 2+1 ))/(L A 2*(6*a A 5*ln(a)*p A 4+6*a A 4*ln(a)*p A 4- 
12*a A 3*ln(a)*p A 4-14*a A 5*p A 4+22*a A 4*p A 4+a A 3*p A 4-13*a A 2*p A 4+5*a*p A 4- 
P A 4+12*a A 3*ln(a)*P A 3-22*a A 3*P A 3+36*a A 2*P A 3-18*a*P A 3+4*P A 3- 
6*a A 4*p A 2+24*a A 3*p A 2-36*a A 2*p A 2+24*a*p A 2-6*p A 2+2*a A 5*p-4*a A 4*p- 
4*a A 3*p+16*a A 2*p-14*a*p+4*p+a A 3-3*a A 2+3*a-1)) 

C=(3*v1 *(a-1 )*a A 3*p*(3*a*p A 2-3*p A 2+2*p- 

a+1 ))/(L*(6*a A 5*ln(a)*p A 4+6*a A 4*ln(a)*p A 4-12*a A 3*ln(a)*p A 4- 
14*a A 5*p A 4+22*a A 4*p A 4+a A 3*p A 4-13*a A 2*p A 4+5*a*p A 4-p A 4+12*a A 3*ln(a)*p A 3- 
22*a A 3*P A 3+36*a A 2*P A 3-18*a*P A 3+4*P A 3-6*a A 4*P A 2+24*a A 3*P A 2- 
36*a A 2*p A 2+24*a*p A 2-6*p A 2+2*a A 5*p-4*a A 4*p-4*a A 3*p+16*a A 2*p- 
14*a*p+4*p+a A 3-3*a A 2+3*a-1)) 

D=-(v1 *(a A 3*p A 4*ln(L A 12*a A 12*p A 12)-a A 3*p A 3*ln(L A 12*a A 12*p A 12)- 
a A 5*p A 4*ln(L A 6*a A 6*p A 6)-a A 4*p A 4*ln(L A 6*a A 6*p A 6)+8*a A 5*p A 4- 
19*a A 4*P A 4+2*a A 3*P A 4+13*a A 2*p A 4-5*a*p A 4+p A 4+16*a A 3*p A 3- 
36*a A 2*P A 3+18*a*P A 3-4*P A 3+3*a A 4*P A 2-21 *a A 3*p A 2+36*a A 2*p A 2- 
24*a*p A 2+6*p A 2-2*a A 5*p+4*a A 4*p+4*a A 3*p-16*a A 2*p+14*a*p-4*p-a A 3+3*a A 2- 
3*a+1))/(6*a A 5*ln(a)*p A 4+6*a A 4*ln(a)*p A 4-12*a A 3*ln(a)*p A 4- 
14*a A 5*p A 4+22*a A 4*p A 4+a A 3*p A 4-13*a A 2*p A 4+5*a*p A 4-p A 4+12*a A 3*ln(a)*p A 3- 
22*a A 3*P A 3+36*a A 2*P A 3-18*a*P A 3+4*P A 3-6*a A 4*P A 2+24*a A 3*P A 2- 
36*a A 2*p A 2+24*a*p A 2-6*p A 2+2*a A 5*p-4*a A 4*p-4*a A 3*p+16*a A 2*p- 
14*a*p+4*p+a A 3-3*a A 2+3*a-1) 

E1=(3*v1*(a-1) A 4*p*(2*p-a-2))/(L*(6*a A 5*ln(a)*p A 4+6*a A 4*ln(a)*p A 4- 

12*a A 3*ln(a)*p A 4-14*a A 5*p A 4+22*a A 4*p A 4+a A 3*p A 4-13*a A 2*p A 4+5*a*p A 4- 
P A 4+12*a A 3*ln(a)*P A 3-22*a A 3*P A 3+36*a A 2*P A 3-18*a*P A 3+4*P A 3- 
6*a A 4*p A 2+24*a A 3*p A 2-36*a A 2*p A 2+24*a*p A 2-6*p A 2+2*a A 5*p-4*a A 4*p- 
4*a A 3*p+16*a A 2*p-14*a*p+4*p+a A 3-3*a A 2+3*a-1)) 

F=((3*v1*a A 4-12*v1*a A 3+18*v1*a A 2-12*v1*a+3*v1)*p A 2+(- 

2*v1 *a A 5+4*v1 *a A 4+4*v1 *a A 3-16*v1 *a A 2+14*v1 *a-4*v1 )*p-v1 *a A 3+3*v1 *a A 2- 
3*v1*a+v1 )/(P A 4*(a A 3*(12*ln(L*(a-1)*p+L*P)-12*ln(L*P)-1)+a A 5*(-6*ln(L*(a- 
1 )*p+L*P)+6*ln(L*P)+14)+a A 4*(-6*ln(L*(a-1 )*p+L*P)+6*ln(L*P)-22)+13*a A 2- 
5*a+1 )+P A 3*(a A 3*(-12*ln(L*(a-1 )*p+L*p)+12*ln(L*p)+22)-36*a A 2+18*a- 
4)+(6*a A 4-24*a A 3+36*a A 2-24*a+6)*P A 2+(-2*a A 5+4*a A 4+4*a A 3-16*a A 2+14*a- 
4)*p-a A 3+3*a A 2-3*a+1) 


De manera que se tiene que: 
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• Para: 0 < x < (3L 


0 2 (x)vi 

AL 3 |3 3 ln(L(3 + x(a — 1)) _ AL 4 (3 4 + (BL 3 - BL 3 a)|3 3 _ 

a 3 — 3a 2 + 3a — 1 + (2La 3 — 6La 2 + 6La — 2L)(3 + x(2a 4 — 8a 3 + 12a 2 — 8a + 2) 

E * Ia 

+ Cx+D (13) 

• Para: f$L <x < L 


Ax 3 Bx 2 

^(»> 1 = 6 Elaa 3 2 +E1^ + F (14) 


Dada la longitud de las funciones y las constantes, el valor de vl puede ser 
simplificado cuando se asignen valores numericos a las variables simbolicas. 

Para el calculo de momento en el nudo inicial, es necesario resolver 
nuevamente el sistema de ecuaciones tomando la condicion de un giro 01 en el 
nudo inicial, mientras los demas desplazamientos son 0. 


x = 0 9 = 61 

x = 0 v = 0 
x—L 9—0 
x — L 17 — 0 


De manera que solo las ecuaciones 7 y 8 cambian mientras el resto 
permanece igual. 


AL 4 B 4 + (BL 3 a - BL 3 )B 3 

r _____ r — qi 

E\a(2L 2 a 2 — 4 L 2 a + 2 L 2 )p 2 

AL 3 P 3 ln (L/t) AL*P 4 + (BL 3 - BL 3 a)/3 3 

a 3 - 3a 2 + 3a - 1 (2 La 3 - 6La 2 + 6 La - 2 L)B 

- k - ~ +D = ° 


Con estas ecuaciones se resuelve el sistema y se obtiene las siguientes 
constantes de integracion. 
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A=-(((6*E*la*a A 7-24*E*la*a A 6+36*E*la*a A 5- 

24*E*la*a A 4+6*E*la*a A 3)*(3 A 2+6*E*la*a A 6-18*E*la*a A 5+18*E*la*a A 4- 
6*E*la*a A 3)*61)/(p A 4*(L A 2*a A 3*(12*ln(L*(a-1)*p+L*p)-12*ln(L*p)-1)+L A 2*a A 5*(- 
6*ln(L*(a-1)*p+L*P)+6*ln(L*P)+14)+L A 2*a A 4*(-6*ln(L*(a-1)*p+L*P)+6*ln(L*P)- 
22)+13*L A 2*a A 2-5*L A 2*a+L A 2)+P A 3*(L A 2*a A 3*(-12*ln(L*(a- 
1)*P+L*P)+12*ln(L*P)+22)-36*L A 2*a A 2+18*L A 2*a-4*L A 2)+(6*L A 2*a A 4- 
24*L A 2*a A 3+36*L A 2*a A 2-24*L A 2*a+6*L A 2)*P A 2+(- 
2*L A 2*a A 5+4*L A 2*a A 4+4*L A 2*a A 3-16*L A 2*a A 2+14*L A 2*a-4*L A 2)*p- 
L A 2*a A 3+3*L A 2*a A 2-3*L A 2*a+L A 2) 


B=-((P A 3*(E*la*a A 6*(-12*ln((L*(a-1)*p+L*P))+12*ln((L*P))+22)- 
36*E*la*a A 5+18*E*la*a A 4-4*E*la*a A 3)-4*E*la*a A 6+12*E*la*a A 5- 
12*E*la*a A 4+4*E*la*a A 3)*01)/(P A 4*(L*a A 3*(12*ln((L*(a-1)*p+L*P))-12*ln((L*P))- 
1 )+L*a A 5*(-6*ln((L*(a-1)*p+L*P))+6*ln((L*P))+14)+L*a A 4*(-6*ln((L*(a- 
1 )*p+L*P))+6*ln((L*P))-22)+13*L*a A 2-5*L*a+L)+p A 3*(L*a A 3*(-12*ln((L*(a- 
1 )*P+L*P))+12*ln((L*P))+22)-36*L*a A 2+18*L*a-4*L)+(6*L*a A 4- 
24*L*a A 3+36*L*a A 2-24*L*a+6*L)*P A 2+(-2*L*a A 5+4*L*a A 4+4*L*a A 3- 
16*L*a A 2+14*L*a-4*L)*p-L*a A 3+3*L*a A 2-3*L*a+L) 

C=-((P A 4*(a A 3*(12*ln((L*(a-1)*p+L*P))-12*ln((L*P))-19)+a A 4*(-18*ln((L*(a- 
1)*P+L*P))+18*ln((L*P))-3)+9*a A 5+18*a A 2-6*a+1)+p A 3*(a A 4*(12*ln((L*(a- 
1)*p+L*p))-12*ln((L*p))-22)+a A 3*(-12*ln((L*(a-1)*p+L*p))+12*ln((L*p))+58)- 
54*a A 2+22*a-4)+(-3*a A 5+24*a A 4-57*a A 3+60*a A 2-30*a+6)*P A 2+(-6*a A 4+22*a A 3- 
30*a A 2+18*a-4)*p+a A 4-4*a A 3+6*a A 2-4*a+1 )*01 )/(p A 4*(a A 4*(18*in((L*(a- 
1 )*P+L*P))-18*ln((L*P))+21)+a A 6*(-6*ln((L*(a-1 )*p+L*p))+6*ln((L*p))+14)+a A 3*(- 
12*ln((L*(a-1)*p+L*P))+12*ln((L*P))+14)-36*a A 5-18*a A 2+6*a- 
1)+P A 3*(a A 3*(12*ln((L*(a-1)*p+L*P))-12*ln((L*P))-58)+a A 4*(-12*ln((L*(a- 
1 )*p+L*P))+12*ln((L*P))+22)+54*a A 2-22*a+4)+(6*a A 5-30*a A 4+60*a A 3- 
60*a A 2+30*a-6)*p A 2+(-2*a A 6+6*a A 5-20*a A 3+30*a A 2-18*a+4)*p-a A 4+4*a A 3- 
6*a A 2+4*a-1) 

D=-((P A 5*(L*a A 6*(6*ln((L*(a-1)*p+L*P))-8)+18*L*a A 4*ln((L*P))+L*a A 3*(-6*ln((L*P))- 

1)+L*a A 5*(9- 

18*ln((L*P))))+p A 3*(L*a A 5*(6*ln((L*P))+3)+L*a A 3*(6*ln((L*P))+3)+L*a A 4*(- 
12*ln((L*P))-6))+(2*L*a A 6-6*L*a A 5+6*L*a A 4- 

2*L*a A 3)*p A 2)*01)/(P A 3*(a A 4*(24*ln((L*(a-1)*p+L*P))-24*ln((L*P))-80)+a A 3*(- 

12*ln((L*(a-1)*p+L*P))+12*ln((L*P))+112)+a A 5*(-12*ln((L*(a- 

1 )*P+L*P))+12*ln((L*P))+22)-76*a A 2+26*a-4)+P A 4*(a A 5*(18*ln((L*(a-1 )*p+L*p))- 

18*ln((L*P))+57)+a A 3*(12*ln((L*(a-1)*p+L*P))-12*ln((L*P))-32)+a A 6*(6*ln((L*(a- 

1)*P+L*P))-6*ln((L*P))-50)+a A 7*(-6*ln((L*(a-1)*p+L*P))+6*ln((L*P))+14)+a A 4*(- 

30*ln((L*(a-1)*p+L*P))+30*ln((L*P))-7)+24*a A 2-7*a+1)+(6*a A 6-36*a A 5+90*a A 4- 

120*a A 3+90*a A 2-36*a+6)*P A 2+(-2*a A 7+8*a A 6-6*a A 5-20*a A 4+50*a A 3- 

48*a A 2+22*a-4)*p-a A 5+5*a A 4-10*a A 3+10*a A 2-5*a+1) 

E1=((P A 3*(a A 3*(-12*ln((L*(a-1)*p+L*P))+12*ln((L*P))+22)-36*a A 2+18*a-4)+(3*a A 4- 
12*a A 3+18*a A 2-12*a+3)*p A 2-a A 3+3*a A 2-3*a+1)*01)/(P A 4*(a A 3*(12*ln((L*(a- 
1 )*P+L*P))-12*ln((L*P))-1 )+a A 5*(-6*ln((L*(a-1 )*p+L*p))+6*ln((L*p))+14)+a A 4*(- 
6*ln((L*(a-1 )*p+L*p))+6*ln((L*p))-22)+13*a A 2-5*a+1 )+p A 3*(a A 3*(-12*ln((L*(a- 
1 )*P+L*P))+12*ln((L*P))+22)-36*a A 2+18*a-4)+(6*a A 4-24*a A 3+36*a A 2- 
24*a+6)*p A 2+(-2*a A 5+4*a A 4+4*a A 3-16*a A 2+14*a-4)*p-a A 3+3*a A 2-3*a+1) 
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F=-((P A 3*(L*a A 3*(-6*ln((L*(a-1)*p+L*P))+6*ln((L*P))+11 )-18*L*a A 2+9*L*a- 

2*L)+(2*L*a A 4-8*L*a A 3+12*L*a A 2-8*L*a+2*L)*p A 2)*61)/(P A 4*(a A 3*(12*ln((L*(a- 
1 )*p+L*P))-12*ln((L*P))-1 )+a A 5*(-6*ln((L*(a-1)*p+L*P))+6*ln((L*P))+14)+a A 4*(- 
6*ln((L*(a-1 )*p+L*p))+6*ln((L*p))-22)+13*a A 2-5*a+1 )+p A 3*(a A 3*(-12*ln((L*(a- 
1 )*P+L*P))+12*ln((L*P))+22)-36*a A 2+18*a-4)+(6*a A 4-24*a A 3+36*a A 2- 
24*a+6)*p A 2+(-2*a A 5+4*a A 4+4*a A 3-16*a A 2+14*a-4)*p-a A 3+3*a A 2-3*a+1) 

Teniendo como resultado las siguientes funciones de forma: 

• Para: 0 < x < fiL 

03 W *®1 

AL 3 [3 3 ln(L(3 + x(« — 1)) _ AL 4 (3 4 + (BL 3 - BL 3 a)|3 3 _ 

a 3 — 3a 2 + 3a — 1 (2La 3 — 6La 2 + 6La — 2L)P + x(2a 4 — 8a 3 + 12a 2 

“ H 

+ Cx + D 

• Para: (?L < x < L 



Ax 3 Bx 2 

03W - 8 ‘ = 4tH + EU + F 06) 

Pasando el analisis al nudo final de la media cartela se tiene que las 
condiciones para determinar la funcion 0 5 (x) son: 

x — 0 9 — 0 

x — 0 v — 0 
x — L 9 — 0 
x — L v — v 2 

Con estas consideraciones, las ecuaciones que cambian para el calculo de 
las constantes de integracion son: 

AL 4 /? 4 + (BL 3 a-BL 3 )/3 3 
C E\a(2L 2 a 2 — 4 L 2 a + 2L 2 )P 2 


AL 3 , BL 2 
6 + 2 
iilaa 3 


+ EIL + F — v 2 
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Resolviendo el sistema lineal se tiene que las constantes son: 

A=((6*E*la*v2*a A 8-12*E*la*v2*a A 7-12*E*la*v2*a A 6+48*E*la*v2*a A 5- 
42*E*la*v2*a A 4+12*E*la*v2*a A 3)*(3+12*E*la*v2*a A 6- 
36*E*la*v2*a A 5+36*E*la*v2*a A 4-12*E*la*v2*a A 3)/((3 A 4*(L A 3*a A 3*(12*ln(L*(a- 
1 )*|3+L*P)-12*ln(L*P)-1 )+L A 3*a A 5*(-6*ln(L*(a- 
1)*P+L*P)+6*ln(L*P)+14)+L A 3*a A 4*(-6*ln(L*(a-1)*p+L*P)+6*ln(L*P)- 
22)+13*L A 3*a A 2-5*L A 3*a+L A 3)+P A 3*(L A 3*a A 3*(-12*ln(L*(a- 
1)*P+L*P)+12*ln(L*P)+22)-36*L A 3*a A 2+18*L A 3*a-4*L A 3)+(6*L A 3*a A 4- 
24*L A 3*a A 3+36*L A 3*a A 2-24*L A 3*a+6*L A 3)*P A 2+(- 
2*L A 3*a A 5+4*L A 3*a A 4+4*L A 3*a A 3-16*L A 3*a A 2+14*L A 3*a-4*L A 3)*p- 
L A 3*a A 3+3*L A 3*a A 2-3*L A 3*a+L A 3) 

B=-((6*E*la*v2*a A 7-24*E*la*v2*a A 6+36*E*la*v2*a A 5- 
24*E*la*v2*a A 4+6*E*la*v2*a A 3)*P A 2+6*E*la*v2*a A 6- 
18*E*la*v2*a A 5+18*E*la*v2*a A 4-6*E*la*v2*a A 3)/(P A 4*(L A 2*a A 3*(12*ln(L*(a- 
1 )*P+L*P)-12*ln(L*P)-1 )+L A 2*a A 5*(-6*ln(L*(a- 
1)*P+L*P)+6*ln(L*P)+14)+L A 2*a A 4*(-6*ln(L*(a-1)*p+L*P)+6*ln(L*P)- 
22)+13*L A 2*a A 2-5*L A 2*a+L A 2)+P A 3*(L A 2*a A 3*(-12*ln(L*(a- 
1)*P+L*P)+12*ln(L*P)+22)-36*L A 2*a A 2+18*L A 2*a-4*L A 2)+(6*L A 2*a A 4- 
24*L A 2*a A 3+36*L A 2*a A 2-24*L A 2*a+6*L A 2)*P A 2+(- 
2*L A 2*a A 5+4*L A 2*a A 4+4*L A 2*a A 3-16*L A 2*a A 2+14*L A 2*a-4*L A 2)*p- 
L A 2*a A 3+3*L A 2*a A 2-3*L A 2*a+L A 2) 


C=((9*v2*a A 5-18*v2*a A 4+9*v2*a A 3)*P A 3+(6*v2*a A 4-6*v2*a A 3)*P A 2+(- 

3*v2*a A 5+6*v2*a A 4-3*v2*a A 3)*P)/(P A 4*(L*a A 3*(12*ln(L*(a-1)*p+L*P)-12*ln(L*P)- 

1)+L*a A 5*(-6*ln(L*(a-1)*p+L*P)+6*ln(L*P)+14)+L*a A 4*(-6*ln(L*(a- 

1 )*p+L*P)+6*ln(L*P)-22)+13*L*a A 2-5*L*a+L)+p A 3*(L*a A 3*(-12*ln(L*(a- 

1)*P+L*P)+12*ln(L*P)+22)-36*L*a A 2+18*L*a-4*L)+(6*L*a A 4-24*L*a A 3+36*L*a A 2- 

24*L*a+6*L)*p A 2+(-2*L*a A 5+4*L*a A 4+4*L*a A 3-16*L*a A 2+14*L*a-4*L)*p- 

L*a A 3+3*L*a A 2-3*L*a+L) 

D=(p A 4*(a A 5*(6*v2*ln(L*p)+6*v2)+a A 4*(6*v2*ln(L*p)-3*v2)+a A 3*(-12*v2*ln(L*p)- 

3*v2))+a A 3*P A 3*(12*v2*ln(L*P)+6*v2)+(3*v2*a A 4- 

3*v2*a A 3)*p A 2)/(P A 4*(a A 3*(12*ln(L*(a-1)*p+L*P)-12*ln(L*P)-1)+a A 5*(-6*ln(L*(a- 
1 )*p+L*P)+6*ln(L*P)+14)+a A 4*(-6*ln(L*(a-1 )*p+L*P)+6*ln(L*P)-22)+13*a A 2- 
5*a+1 )+p A 3*(a A 3*(-12*ln(L*(a-1 )*p+L*p)+12*ln(L*p)+22)-36*a A 2+18*a- 
4)+(6*a A 4-24*a A 3+36*a A 2-24*a+6)*P A 2+(-2*a A 5+4*a A 4+4*a A 3-16*a A 2+14*a- 
4)*p-a A 3+3*a A 2-3*a+1) 

E1 =((6*v2*a A 4-24*v2*a A 3+36*v2*a A 2-24*v2*a+6*v2)*P A 2+(- 
3*v2*a A 5+6*v2*a A 4+6*v2*a A 3-24*v2*a A 2+21*v2*a- 
6*v2)*p)/(p A 4*(L*a A 3*(12*ln(L*(a-1)*p+L*p)-12*ln(L*p)-1)+L*a A 5*(-6*ln(L*(a- 
1)*P+L*P)+6*ln(L*P)+14)+L*a A 4*(-6*ln(L*(a-1)*p+L*P)+6*ln(L*P)-22)+13*L*a A 2- 
5*L*a+L)+p A 3*(L*a A 3*(-12*ln(L*(a-1)*p+L*p)+12*ln(L*p)+22)-36*L*a A 2+18*L*a- 
4*L)+(6*L*a A 4-24*L*a A 3+36*L*a A 2-24*L*a+6*L)*P A 2+(- 
2*L*a A 5+4*L*a A 4+4*L*a A 3-16*L*a A 2+14*L*a-4*L)*p-L*a A 3+3*L*a A 2-3*L*a+L) 

F=(P A 4*(a A 3*(v2*(12*ln(L*(a-1)*p+L*P)-1)-12*v2*ln(L*P))+a A 5*(v2*(14-6*ln(L*(a- 
1 )*p+L*P))+6*v2*ln(L*P))+a A 4*(v2*(-6*ln(L*(a-1 )*p+L*p)- 
22)+6*v2*ln(L*p))+13*v2*a A 2-5*v2*a+v2)+p A 3*(a A 3*(v2*(22-12*ln(L*(a- 
1 )*P+L*P))+12*v2*ln(L*P))-36*v2*a A 2+18*v2*a-4*v2)+(3*v2*a A 4- 
12*v2*a A 3+18*v2*a A 2-12*v2*a+3*v2)*p A 2)/(P A 4*(a A 3*(12*ln(L*(a-1)*p+L*P)- 
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12*ln(L*P)-1)+a A 5*(-6*ln(L*(a-1)*p+L*P)+6*ln(L*P)+14)+a A 4*(-6*ln(L*(a- 
1 )*p+L*P)+6*ln(L*P)-22)+13*a A 2-5*a+1 )+p A 3*(a A 3*(-12*ln(L*(a- 
1)*P+L*P)+12*ln(L*P)+22)-36*a A 2+18*a-4)+(6*a A 4-24*a A 3+36*a A 2- 
24*a+6)*p A 2+(-2*a A 5+4*a A 4+4*a A 3-16*a A 2+14*a-4)*p-a A 3+3*a A 2-3*a+1) 

Con estas constantes de integracion las funciones resultantes en funcion de v 2 

son: 


• Para: 0 < x < (SL 


05 (X)V 2 

AL 3 (S 3 ln(L(3 + x(« — 1)) _ AL 4 (3 4 + (BL 3 - BL 3 a)p 3 _ 

a 3 — 3a 2 + 3a — 1 (2La 3 — 6La 2 + 6La — 2L)P + x(2a 4 — 8a 3 + 12a 2 — 8a + 2) 

E * Ia 


+ Cx + D 


(17) 


• Para: [SL < x < L 


Ax 3 Bx 2 

V s (.x)v 2 = 6 Elaa3 2 +E1 x + F (18) 

Para finalizar el analisis es necesario contemplar las condiciones de la 
funcion 0 6 (x) con las siguientes condiciones: 


o 

II 

o 

II 

<£> 

o 

II 

* 

o 

II 

x = L 

9 = 9 2 

II 

* 

o 

II 


De tal manera que las ecuaciones que cambian son: 


4L 2 


+ BL 


E I n a 3 


■+ E1 = 02 


AL 3 ^ BL 2 

- 4 — - 1 —h E1 L + F = 0 

E I a a 3 


Con las mismas ecuaciones que se ha ido trabajando, se resuelve el 
sistema lineal y se tiene que las constantes de integracion son: 
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A=(P*(-6*E*la*a A 8*02+12*E*la*a A 7*02+12*E*la*a A 6*02- 

48*E*la*a A 5*02+42*E*la*a A 4*02-12*E*la*a A 3*02)+P A 2*(6*E*la*a A 7*02- 
24*E*la*a A 6*02+36*E*la*a A 5*02-24*E*la*a A 4*02+6*E*la*a A 3*02)- 
6*E*la*a A 6*02+18*E*la*a A 5*02- 

18*E*la*a A 4*02+6*E*la*a A 3*02)/(P A 4*(L A 2*a A 3*(12*ln(L*(a-1 )*p+L*p)- 

12*ln(L*P)-1)+L A 2*a A 5*(-6*ln(L*(a-1)*p+L*P)+6*ln(L*P)+14)+L A 2*a A 4*(-6*ln(L*(a- 

1 )*p+L*P)+6*ln(L*P)-22)+13*L A 2*a A 2-5*L A 2*a+L A 2)+p A 3*(L A 2*a A 3*(-12*ln(L*(a- 

1)*P+L*P)+12*ln(L*P)+22)-36*L A 2*a A 2+18*L A 2*a-4*L A 2)+(6*L A 2*a A 4- 

24*L A 2*a A 3+36*L A 2*a A 2-24*L A 2*a+6*L A 2)*P A 2+(- 

2*L A 2*a A 5+4*L A 2*a A 4+4*L A 2*a A 3-16*L A 2*a A 2+14*L A 2*a-4*L A 2)*p- 

L A 2*a A 3+3*L A 2*a A 2-3*L A 2*a+L A 2) 


B=(p A 3*(E*la*a A 6*(-12*ln(L*(a-1)*p+L*p)*02+12*in(L*p)*02+22*02)- 

36*E*la*a A 5*02+18*E*la*a A 4*02-4*E*la*a A 3*02)+P A 2*(6*E*la*a A 7*02- 

24*E*la*a A 6*02+36*E*la*a A 5*02- 

24*E*la*a A 4*02+6*E*la*a A 3*02)+2*E*la*a A 6*02-6*E*la*a A 5*02+6*E*la*a A 4*02- 
2*E*la*a A 3*02)/(P A 4*(L*a A 3*(12*ln(L*(a-1 )*p+L*p)-12*ln(L*p)-1 )+L*a A 5*(- 
6*ln(L*(a-1)*p+L*P)+6*ln(L*P)+14)+L*a A 4*(-6*ln(L*(a-1)*p+L*P)+6*ln(L*P)- 
22)+13*L*a A 2-5*L*a+L)+P A 3*(L*a A 3*(-12*ln(L*(a-1 )*p+L*p)+12*ln(L*p)+22)- 
36*L*a A 2+18*L*a-4*L)+(6*L*a A 4-24*L*a A 3+36*L*a A 2-24*L*a+6*L)*P A 2+(- 
2*L*a A 5+4*L*a A 4+4*L*a A 3-16*L*a A 2+14*L*a-4*L)*p-L*a A 3+3*L*a A 2-3*L*a+L) 

C=(P A 4*(a A 6*(-6*ln(L*(a-1)*p+L*P)*02+6*ln(L*P)*02+14*02)-27*a A 5*02+18*a A 4*02- 

5*a A 3*02)+p*(a A 6*02-3*a A 5*02+3*a A 4*02-a A 3*02)+P A 3*(-9*a A 6*02+27*a A 5*02- 

27*a A 4*02+9*a A 3*02)+P A 2*(-3*a A 5*02+6*a A 4*02- 

3*a A 3*02))/(P A 4*(a A 4*(18*ln(L*(a-1 )*p+L*p)-18*ln(L*p)+21 )+a A 6*(-6*ln(L*(a- 
1 )*p+L*P)+6*ln(L*P)+14)+a A 3*(-12*ln(L*(a-1 )*p+L*p)+12*ln(L*p)+14)-36*a A 5- 
18*a A 2+6*a-1 )+p A 3*(a A 3*(12*ln(L*(a-1 )*p+L*p)-12*in(L*p)-58)+a A 4*(-12*ln(L*(a- 
1)*P+L*P)+12*ln(L*P)+22)+54*a A 2-22*a+4)+(6*a A 5-30*a A 4+60*a A 3- 
60*a A 2+30*a-6)*p A 2+(-2*a A 6+6*a A 5-20*a A 3+30*a A 2-18*a+4)*p-a A 4+4*a A 3- 
6*a A 2+4*a-1) 

D=(p A 5*(L*a A 6*(6*ln(L*(a-1)*p+L*p)*02-8*02)+L*a A 3*(-6*ln(L*p)*02- 

02)+L*a A 5*(9*02- 

18*ln(L*P)*02)+18*L*a A 4*ln(L*P)*02)+p A 4*(L*a A 5*(18*ln(L*P)*02- 

9*02)+L*a A 3*(12*ln(L*P)*02+3*02)+L*a A 6*(6*ln(L*P)*02+15*02)+L*a A 7*(- 

6*ln(L*P)*02-6*02)+L*a A 4*(-3O*ln(L*P)*02- 

3*02))+P A 3*(L*a A 4*(12*ln(L*P)*02+6*02)+L*a A 5*(-6*ln(L*P)*02-3*02)+L*a A 3*(- 

6*ln(L*P)*02-3*02))+P A 2*(-L*a A 6*02+3*L*a A 5*02- 

3*L*a A 4*02+L*a A 3*62))/(P A 3*(a A 4*(24*ln(L*(a-1)*p+L*P)-24*ln(L*P)-8O)+a A 3*(- 
12*ln(L*(a-1 )*p+L*p)+12*ln(L*p)+112)+a A 5*(-12*ln(L*(a- 
1)*P+L*P)+12*ln(L*P)+22)-76*a A 2+26*a-4)+p A 4*(a A 5*(18*ln(L*(a-1)*p+L*P)- 
18*ln(L*P)+57)+a A 3*(12*ln(L*(a-1 )*p+L*p)-12*ln(L*p)-32)+a A 6*(6*ln(L*(a- 
1)*P+L*P)-6*ln(L*P)-50)+a A 7*(-6*ln(L*(a-1)*p+L*P)+6*ln(L*P)+14)+a A 4*(- 
30*ln(L*(a-1)*p+L*P)+30*ln(L*P)-7)+24*a A 2-7*a+1)+(6*a A 6-36*a A 5+90*a A 4- 
120*a A 3+90*a A 2-36*a+6)*P A 2+(-2*a A 7+8*a A 6-6*a A 5-20*a A 4+50*a A 3- 
48*a A 2+22*a-4)*p-a A 5+5*a A 4-10*a A 3+10*a A 2-5*a+1) 
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E1=(P A 4*(a A 3*(12*ln(L*(a-1)*p+L*P)*02-12*ln(L*P)*02-02)+a A 5*(-6*ln(L*(a- 

1 )*p+L*P)*02+6*ln(L*P)*02+14*02)+a A 4*(-6*ln(L*(a-1)*p+L*P)*02+6*ln(L*P)*02- 
22*02)+13*a A 2*02-5*a*02+02)+p*(a A 5*02-2*a A 4*02-2*a A 3*02+8*a A 2*02- 
7*a*02+2*02)+P A 2*(-3*a A 4*02+12*a A 3*02-18*a A 2*02+12*a*02- 
3*02))/(P A 4*(a A 3*(12*ln(L*(a-1)*p+L*P)-12*ln(L*P)-1)+a A 5*(-6*ln(L*(a- 
1 )*p+L*P)+6*ln(L*P)+14)+a A 4*(-6*ln(L*(a-1 )*p+L*P)+6*ln(L*P)-22)+13*a A 2- 
5*a+1 )+p A 3*(a A 3*(-12*ln(L*(a-1 )*p+L*p)+12*ln(L*p)+22)-36*a A 2+18*a- 
4)+(6*a A 4-24*a A 3+36*a A 2-24*a+6)*P A 2+(-2*a A 5+4*a A 4+4*a A 3-16*a A 2+14*a- 
4)*p-a A 3+3*a A 2-3*a+1) 

F=-(P A 4*(L*a A 3*(12*ln(L*(a-1 )*p+L*p)*02-12*ln(L*p)*02-02)+L*a A 5*(-6*ln(L*(a- 
1 )*p+L*P)*02+6*ln(L*P)*02+14*02)+L*a A 4*(-6*ln(L*(a- 

1)*P+L*P)*02+6*ln(L*P)*02-22*02)+13*L*a A 2*02-5*L*a*02+L*02)+P A 3*(L*a A 3*(- 
6*ln(L*(a-1 )*p+L*p)*02+6*ln(L*p)*02+11 *02)-18*L*a A 2*02+9*L*a*02- 
2*L*02)+P A 2*(L*a A 4*02-4*L*a A 3*02+6*L*a A 2*02- 

4*L*a*02+L*02))/(P A 4*(a A 3*(12*ln(L*(a-1 )*p+L*p)-12*ln(L*p)-1 )+a A 5*(-6*ln(L*(a- 
1 )*p+L*P)+6*ln(L*P)+14)+a A 4*(-6*ln(L*(a-1 )*p+L*P)+6*ln(L*P)-22)+13*a A 2- 
5*a+1 )+P A 3*(a A 3*(-12*ln(L*(a-1 )*p+L*p)+12*ln(L*p)+22)-36*a A 2+18*a- 
4)+(6*a A 4-24*a A 3+36*a A 2-24*a+6)*P A 2+(-2*a A 5+4*a A 4+4*a A 3-16*a A 2+14*a- 
4)*p-a A 3+3*a A 2-3*a+1) 

De la misma manera las funciones de forma resultantes en funcion de 0 2 son: 

• Para: 0 < x < pL 


0 6 (x) * 0 2 

AL 3 |3 3 ln(L(3 + x(a — 1)) _ AL 4 (3 4 + (BL 3 - BL 3 a)|3 3 _ 

a 3 — 3a 2 + 3a — 1 (2La 3 — 6La 2 + 6La — 2L)(3 + x(2a 4 — 8a 3 + 12a 2 — 8a + 2) 

+ Cx + D 


e*l 


(19) 




Para: (iL < x < L 


0 6 (%) * e 2 


Ax 3 ^ Bx 2 

_6_2_+ pivXC 

El a a 3 


( 20 ) 


Cabe destacar que, por la longitud de las constantes, no se ha hecho la 
simplificacion del modulo de elasticidad multiplicado por la inercia I a , por lo que, al 
momento de realizar el calculo manual, los valores que estas tomen son 
despreciables. 


3. TABLAS DE CALCULO 

Con las funciones de forma 13 a 20, se han generado tablas con las relaciones 
a y p mas comunes para encontrar las reacciones de empotramiento perfecto 
considerando una carga uniformemente distribuida a lo largo de toda la viga. 
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Tablas de cortante 

El coeficiente presentado en las tablas se da para ser usado de la siguiente 
manera: 


V = C*P*L 

Siendo, C el coeficiente determinado en la tabla, P el valor de carga aplicado y 
L la longitud de la viga. 

Tabla 1 Coeficientes para determinar cortante en el nudo inicial de la viga. 


p.a 

0,05 

0.9 

0.95 

0,8 

0,75 

mm 


... 

0.65 

0,5 

0,45 0.4 

mn 

a.50519 

0.51078 

0,51684 

0.52343 

0.53061 

0.53848 

0.54716 

0,55676 

0.56743 

0,57933 

0,59264 Q.607S4 

0.45 

0.50522 

0 51006 

0,51695 

0.5235$ 

0.53080 

0 5386$ 

0.54733 

0.55683 

0,56729 

0.57082 

0.59152 0.60540 

0,4 

0,50524 

0.51084 

0,51697 

0.52356 | 0.5JD/1 


0,54691 

0,55609 

0,56608 


0.58866 0.601JJ 

0,35 

0,50520 

0,51078 

0,51678 

0.52322 

0.53016 


0.54565 

055428 

0,56333 


0.58392 0 r 59498 

E3 

0.505« 

0,51050 

0.51627 

0.52242 

0.5J597 

0.53594 

0.54333 

055117 

0.55943 




0.50444 

0,50994 

0,51533 

0.52100 

0,52696 

0,53322 

0,53971. 

0,54651 

0,55 564 


nmnm 

mn 

0,50442 

0.50903 

0,51383 

0.51881 

0,52397 

052930 

0.53479 

0,34041 

0,54613 


IIIIIIWII 

0,1 s 

0.50375 

0,50767 

0.51165 

0.51572 

0.51987 

0.5240$ 

0.52934 

053262 

0,53691 

0,5411$ 

0.54540 0.54954 | 


Tabla 2 Coeficientes para determinar el cortante en el nudo final de la viga. 


(Sia 

0.96 

0.9 

0.85 | 

0.8 

0.75 

mm 

0.55 

0.6 

0.56 

mm 

0.45 

0.4 

0.5 

0.49401 

0.40922 





0.45284 

0.44324 

043257 

0.42067 

0.40736 

0.39246 

0,45 

0.49478 

0,48914 

0.48304 

0,47642 

0,46920 

0/46132 




0,42118 

0.40B48 

0.394S2 

0.4 

0.49476 

0,48912 

0.4S303 

0,47644 

0.46929 

0,46153 

0,45309 

0,44391 

0,43392 



0,19867 

0,26 

0.4948C 

0.48922 

0,48322 

0,47678 

□.46984 

0.46238 

0,45435 

0,44572 

0,43647 


041608 

0,40502 

0,3 

0.49492 

0,48950 

0,48373 

0,47758 

0.47103 

0,46406 

0,45667 

0,44883 

0,44057 


042286 

0.413S3 

0.25 

0.49516 

0.49006 





0.46024 

0,45342 

0,44636 



0,42408 

02 

0.49558 

0.49097 

0,48617 

0,48119 

0,47603 

0,47070 

0,46521 

0,45959 

0,45387 

C.4480S 

0.44225 

0,43646 

0,16 

0.49622 

0.49233 

0.4*835 | 

0.48428 

0.40013 

0/47592 

0.47166 

0.46738 

046309 

0.45882 

045460 

0.45046 


En caso de que la cartela se encuentre en el nudo final de la viga, los 
coeficientes para determinar cortante en el nudo inicial deberan ser tomados de la 
tabla 2 y para el cortante en el nudo final se debe tomar los coeficientes de la tabla 
1 . 


Tablas para momento 

Los coeficientes presentados en estas tablas se deben utilizar de la siguiente 
manera: 


M — C * P * L 2 

Siendo, C el coeficiente determinado en la tabla, P el valor de carga aplicado y 
L la longitud de la viga. 
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Tabla 3 Coeficientes para determinar el momento en el nudo inicial. 



0.95 

0.9 

086 

08 

0,75 

0.7 

0.66 

0,6 


0.6 

0.46 

wm 

0.5 

0,08676 

0,09<M9 

0.094S9 

0.09910 

0,10*10 

0.10960 

0.11570 

0,17760 

0.13030 

0,13900 

a.14880 

0.15990 

0.46 

0,08685 







0,17790 

0.13050 

0,13890 


0,15850 

0.4 

0.08690 

0,09077 

0.09500 

0.09960 

0,10460 

0 11000 

0.11600 

0.17760 

0.17980 

0.13760 

0,14670 

0,15540 

0.36 

0.08690 

0.09075. 

O.OM'W 

0 09940 

0.10420 

0.10950 

0.11520 

0.17130 

0.17790 

0.13500 

0.14760 

0.15060 

0,1 

0.08683 

0.09057 

0,0946(1 

009880 

0.10340 

0,10830 

0.11350 

011900 

0.17490 

0.13110 

ai3750 

0.14420 

0 « 

0.08666 

0,09018 

0,09390 

0.09780 

0.10200 

010630 

0.11090 

0.11570 

0,12060 

0.17580 

0.13100 

0.13640 

mn 

0.08636 

0.08953 

0.092*1 

009630 

0.09090 

0.10350 

0.10730 

0.11120 

0.11520 

0.11920 

0.12330 

0.12740 

0.16 

0.08591 

0.08457 

0,091X1 

0J09410 

0.09690 

0J&9980 

0.10270 

a10570 

0.10860 

0.11160 

0.11450 

0.11740 


Tabla 4 Coeficientes para determinar el momento en el nudo final. 




Bi 

• 99 

M 

8.79 

Bi 

0*9 

13 

0 64 

0.S 

o,*» 

M 

08 

008157 


0,07775 

0.07547 

0,07345 

007109 

0,06855 

0.04581 

0,04284 

0JOM65 

0.05817 

005238 

0.48 

008162 

-007181 

01077*9 

007584 | 

007346 

4)07132 

OOM82 

04M12 

006321 

4)1)8007 

0 05689 

-005306 

04 

-0.Q81M 

wrnro 

0,07800 

007600 | 

007310 

4)07156 

006911 1 

000649 

O04844 

4)06064 

-0 097 M 

4)05410 

024 

008170 


a07813 

0.07617 

007409 

0/1*7185 

0.04953 

0.04704 

0.04439 

00614) 


005364 

0.3 

006175 

0/38007 

O0761U 

007643 

007445 

4)/37237 

0,0701* 

4>Ab78* 

006647 

4)08297 

4106039 

-006775 

0 29 

-0 06187 

■amm 

007117 

0076» ! 

007901 

0/27112 

007115 ! 

000911 

0067011 

41064*3 

-006267 

00804* 

0.2 

0 0615*4 

Q/XK&O 

O.079X) 

0,07745 1 

007585 

007421 

0,07253 ' 0,07081 

0,06907 

006732 

0JO6557 

006383 

8.19 

0062 U 

0j08080 

O07964 

087635 

007704 

007572 

007439 0.07306 

007172 

4)07060 

4L08905 

-0/36782 


En caso de que la cartela se encuentre en el nudo final de la viga, los 
coeficientes para determinar el momento en el nudo inicial deberan ser tomados de 
la tabla 4 y para el momento en el nudo final se debe tomar los coeficientes de la 
tabla 3, conjuntamente se debe invertir los signos en el calculo. 

4. NUEVAS FUNCIONES DE CEINCI-LAB 

En Aguiar y Bonilla (2020) se presentaron las funciones para el calculo de 
la matriz de rigidez, de una viga con cartelas, para tres condiciones a saber: i) Viga 
con cartela sin nudo intermedio; ii) Cartela al lado izquierdo; iii) Cartela al lado 
derecho. 

En este artfculo se ha presentado el calculo de las acciones de 
empotramiento perfecto para dos casos: i) Cartela al lado izquierdo; ii) Cartela al 
lado derecho, que se muestran en la figura 4 b. 

Las tablas 1 a 4, sirven para calcular las acciones de empotramiento 
perfecto para una cartela a la izquierda o una cartela a la derecha, ante una carga 
uniforme distribuida. Ahora, en lugar de usar las tablas se tiene el siguiente 
programa. 
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• [REP]=EP_media_cartela(b,h,H,L,Lv,carga,Tipo) 

El programa EP_media_cartela, encuentra las acciones de empotramiento 
perfecto para una viga con cartela a la izquierda o a la derecha; el codigo para saber 
la ubicacion de la cartela es mediante la variable Tipo, si la cartela esta a la izquierda 
Tipo=0, si esta a la derecha Tipo=1. 

Los restantes datos son la base b\ la altura en el centro de luz h\ la altura 
en el extremo de la viga H\ la longitud de la viga L\ la longitud de la cartela L v \ el 
valor de la carga uniforme distribuida se da a continuacion en la variable carga, 
finalmente se indica la ubicacion de la cartela en Tipo. 

Las acciones de empotramiento perfecto vienen en el vector REP 

• [FF]=fuerzas_cartelas(ngl,ELEM,L,seno,coseno, VC,E,q,Q2) 

Este programa determina las fuerzas y momentos en los elementos con 
cartela al lado izquierda o al lado derecho. Los datos de entrada, son ngl, numero 
de grados de libertad; ELEM, matriz que contiene para cada elemento con cartela: 
b h H Lv Tipo, que son la base, altura en centro de luz, altura en el extremo, la 
longitud de la cartela y la variable Tipo, ya explicada en la funcion anterior; L es un 
vector con la longitud de cada elemento con cartela; seno,coseno, son vectores con 
los valores de las funciones trigonometricas seno y coseno del elemento; VC, es una 
matriz con los vectores de colocacion de los elementos; E, es el modulo de 
elasticidad del hormigon; q, es un vector con los desplazamientos y nudos de la 
estructura; Q 2 , es una matriz con las acciones de empotramiento perfecto de las 
vigas con cartelas. 

5. USO DE FUNCIONES DE CEINCI-LAB 

Se presenta un programa, utilizando el sistema de computacion CEINCI- 
LAB para encontrar las fuerzas y momentos en los elementos de la estructura de 
hormigon armado indicada en la figura 5 a, las cartelas son de 40/70 cm, en el centro 
de luz; 40/100 cm, en los extremos y la longitud de las cartelas lineales simetricas 
es de 200 cm. Las columnas son de 60/120 cm. El modulo E = 2400000 T/m 2 . 

En la figura 5 b, se presenta la numeracion de nudos y elementos, primero 
se numeran las columnas y luego las vigas con cartelas. La matriz de rigidez de la 
estructura se obtiene para los grados de libertad mostrados en la figura 5 c, y se 
encuentra primero la contribucion de las columnas (elementosl a la 9) con la funcion 
krigidez, despues se halla la contribucion a la matriz de rigidez de las cartelas con 
la funcion krigidez_cartela_nudointermedio (elementos 10 al 21) Se recomienda 
ver Aguiar y Bonilla (2020), Aguiar et al. (2019), donde se halla explicado el calculo 
de la matriz de rigidez. 

El marco teorico del Vector de Cargas Q, de la matriz de rigidez de la 
estructura K, del vector de coordenadas q, y de las fuerzas y momentos en cada 
uno de los elementos, se encuentra bien detallado en Aguiar (2020). Esta el marco 
teorico y el uso de las funciones de CEINCI-LAB para los temas que se tratan en el 
libro. 
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Figura 5 Portico plano, con vigas que tienen cartelas lineales simetricas; a) 
Dimensiones de los elementos y de la estructura; b) Numeracion de 

nudos y elementos; c) Identificacion de los grados de libertad. 
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Los programas del sistema de computacion CEINCI-LAB que se utilizan 
para realizar el analisis estatico y dinamico de estructuras con cartelas lineales, se 
hallan en la particion del primero de los autores en research gate, en la carpeta de 
tablas. 


En el Anexo 1, se lista el programa utilizado para resolver la estructura de 
la figura 5 a, y en las tablas 5 y 6 se indican las fuerzas y momentos en coordenadas 
locales de las columnas y de las vigas, en el siguiente orden: fuerza axial en nudo 
inicial N1 en Ton; fuerza de corte en nudo inicial VI en Ton; momento en nudo inicial 
M1 en Ton m; fuerza axial en nudo final N2 en Ton; fuerza de corte en nudo final V2 
en Ton; momento en nudo final M2 en Ton m. 


Tabla 5 Fuerzas y Momentos en coordenadas locales en columnas. 


TIPO 

NUMERO 

N1 

VI 

M1 

N2 

V2 

M2 

COLUMNA 

1 

42.33 

3.41 

34.19 

-42.33 

-3.41 

-21.90 

2 

92.32 

6.03 

37.37 

-92.32 

-6.03 

-15.66 

3 

45.35 

11.55 

43.67 

-45.35 

-11.55 

-2.08 

4 

27.89 

-1.78 

2.23 

-27.89 

1.78 

-8.63 

5 

61.72 

2.36 

9.00 

-61.72 

-2.36 

-0.49 

6 

30.39 

9.42 

22.02 

-30.39 

-9.42 

11.88 

7 

13.86 

-7.06 

-8.90 

-13.86 

7.06 

-16.53 

8 

31.01 

-1.42 

-2.33 

-31.01 

1.42 

-2.80 

9 

15.12 

8.49 

9.89 

-15.12 

-8.49 

20.67 


Tabla 6 Fuerzas y momentos en coordenadas locales en vigas. 


TIPO 

NUMERO 

N1 

VI 

M1 

N2 

V2 

M2 

VIGA 

10 

-1.19 

14.44 

19.66 

1.19 

0.56 

15.04 

11 

-1.19 

-0.56 

-15.04 

1.19 

15.56 

-19.49 

12 

-4.86 

15.04 

26.15 

4.86 

-0.04 

18.49 

13 

2.14 

0.04 

-18.50 

-2.14 

14.96 

-19.94 

14 

-5.29 

14.02 

17.54 

5.29 

0.98 

15.08 

15 

4.71 

-0.98 

-15.07 

-4.71 

15.98 

-21.53 

16 

0.93 

14.73 

24.35 

-0.93 

0.27 

18.48 

17 

0.93 

-0.27 

-18.48 

-0.93 

15.27 

-21.77 

18 

7.06 

13.86 

16.53 

-7.06 

1.14 

15.28 

19 

7.06 

-1.14 

-15.28 

-7.06 

16.14 

-22.14 

20 

8.49 

14.87 

24.93 

-8.49 

0.13 

18.74 

21 

8.49 

-0.13 

-18-74 

-8.49 

15.13 

-20.67 


6. VALIDACION DE RESULTADOS 

Para validar el calculo de las acciones de empotramiento perfecto, en vigas 
con cartela a la izquierda o la derecha, se analiza una viga de 10 m, de longitud; 
cuyas dimensiones son: b = 0.4 m; h = 0.6 m\ H = 0.8 m-, L v = 2.0 m; L = 5.0 m; la 
longitud es 5 m, ya que el modelo de calculo considera nudo en la mitad de la viga. 
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Se calcula de tres maneras, a saber: i) Mediante el uso de las tablas 1 a 4; ii) Con 
el programa Sap 2000; iii) Con el programa EP_media_cartela 

• Uso de tablas 1 a 4 presentadas en este articulo 

2 0.6 
^-= 0.4 « = ^= 0.75 

De las tablas (1 a 4) se obtiene que los coeficientes necesarios para calcular 
las acciones de empotramiento perfecto, estos son: 

Nudo inicial 


CVl = 0.53071 CMl = 0.1046 


Nudo final 


CV2 = 0.46929 CM 2 = -0.07385 

Con estos coeficientes se calculan las reacciones de empotramiento 
perfecto. 

Nudo inicial 


Cortante — 0.53071 * 8 * 5 = 21.23 Tn 
Momento — 0.1046 * 8 * 5 2 = 20.92 Tn*m 


Nudo final 


Cortante — 0.46929 * 8 * 5 = 18.77 Tn 
Momento — —0.07385 * 8 * 5 2 = —14.77 Tn*m 

Para considerar la cartela en el nudo final se toman las consideraciones 
dadas en las tablas. 

• Uso de Sap 2000 

Para el Sap 2000, la viga con cartela a la izquierda se considera como que 
fueran dos vigas, una con variacion de altura lineal y la otra de seccion constante, 
en la figura 6 a, se presentan estos dos elementos y en la figura 6.b, a la izquierda 
se muestra la entrada de datos de la viga con variacion de altura lineal y a la derecha 
los datos de la viga de seccion constante. 

En la figura 6 c, se presenta el diagrama de corte; se ve que los valores son 
21.25 T, en el nudo inicial y 18.75 en el centro de luz. Con el uso de las tablas se 
hallo 21.23 T, en el nudo inicial y 18.77 T., en el nudo final. 

Finalmente, en la figura 6 d, se muestra el diagrama de momentos, hallados 
con Sap 2000, sus valores son 21 T m, en nudo inicial y 14.75 T m, en nudo final. 
Con el uso de tablas se hallo 20.92 T m, en nudo inicial y 14.77 T m, en nudo final. 
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II M II 

A I ( B I C 




(c) 





I 



(d) 

Figura 6 Analisis con Sap 2000; a) Viga conformada por dos elementos; b) Ingreso 
de datos de seccion de elemento con altura variable a la izquierda y de 
altura constante a la derecha; c) Diagrama de Corte que reporta Sap 2000; 
d) Diagrama de Momentos que se halla con Sap 2000. 

• Con programa EP_media_carte!a 

En el apartado 4, se indico el uso del programa EP_media_cartela con esta 
acotacion a continuacion se presenta la entrada de datos y los resultados que arroja 
el programa. 
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[REP]=EP_media_cartela(base,h,H,L,Lv,carga,tipo) 

[REP]=EP_media_cartela(0.4,0.6,0.8,5,2,8,0) 

REP = [21.2282; 20.9117; 18.7718; -14.7705] 

Que son el cortante en nudo inicial, momento en nudo inicial, cortante en 
nudo final y momento en nudo final. 

Las acciones de empotramiento perfecto, halladas con el uso de las tablas 
1 a 4, con Sap 2000 y con la funcion EP_media_cartela son practicamente iguales. 


7. COMENTARIOS Y CONCLUSIONES 

Se han obtenido funciones de forma, para encontrar acciones de 
empotramiento perfecto, ante cualquier tipo de carga, de vigas con cartela a la 
izquierda o a la derecha, ya que se considero nudo en la mitad de la viga. 
Posteriormente se presentaron tablas para encontrar las fuerzas y momentos de 
empotramiento perfecto, para el caso de carga uniforme distribuida. 

Finalmente, se indican nuevas funciones para el sistema de computacion 
CEINCI-LAB que sirven para calcular las cargas de empotramiento perfecto en 
vigas con cartelas ante carga uniforme distribuida y se ha visto su uso mediante el 
calculo de un portico plano con vigas con cartelas lineales simetricas. 

Es una muy buena opcion para cuando se tienen grandes luces disenar la 
estructura con vigas con cartelas y se aspira con este articulo haber aportado a 
entender como se resuelven estas estructuras ante cargas estaticas y lo mas 
importante que el lector conozca la teoria que conduce al calculo de las funciones 
de forma para vigas con cartelas. 
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ANEXO 1 

% Ej_REP_Cartelas 
clear all; clc 

sv=[10.0;12.0]; % luz del vano en metros 
sp=[3.6;3.6;3.6]; % altura de cada piso en metros. 
nod=18;nr=3;E=2400000; 

%% Datos generales 

[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp); 
[X,Y]=glinea_portico2(nv,np,sv,sp,nod,nr); %coordenadas 
[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg); 

%% Geometria de la estructura 
[L,seno,coseno]=longitud (X,Y,NI,NJ); 
for i=1 :nr 

RES(i,:)=[i 1 1 1]; 
end 

[CG,ngl]=cg (nod,nr,RES); 

%dibujogdl(X,Y,NI,NJ,CG); 

%dibujoplano(X,Y,NI,NJ) 

[VC]=vc(NI,NJ,CG); 

% Contribucion de columnas rectangulares 

LC=L(1:9);senC=seno(1:9);cosC=coseno(1:9);VCC=VC(1:9,:); 

Seccion=[1 0.6 1.2 8 1]; [ELEMC]=gelem_portico(Seccion); % tiene base y altura 
[KC]=krigidez(ngl,ELEMC,LC,senC,cosC,VCC,E); 

% Contribucion de vigas con cartelas 

b1=0.40;h1=0.70;H1 =1,0;LV1 =2.0;% Son datos de viga con cartela 
SECCION=[b1 hl H1 LV1]; 

UBICACION=[1 1024 0; 111241; 

1 12 2 4 0; 1 13 2 4 1]; 

[ELEMVC]=gelem_portico_cartela(SECCION,UBICACION); 

LVC=L(10:21 );senVC=seno(10:21 );cosVC=coseno(10:21 );VCVC=VC(10:21,:); 

[KVC]=krigidez_cartela_nudointermedio(ngl,ELEMVC,LVC,senVC,cosVC,VCVC,E) 

KT=KC+KVC; 

%% Fuerzas de los elementos 
njc=3; F=[4 4 0 0; 7 7 0 0;10 10 0 0]; 
nmc=12; Fm=[10 3 4 5 2; 11 3452]; 

[REP1 ]=EP_media_cartela(b1 ,h1 ,H1 ,L(10),LV1 ,Fm(1,2),0);% Empotramientos 
[REP2]=EP_media_cartela(b1 ,h1 ,H1 ,L(11 ),LV1 ,Fm(2,2),1); 
datosp1=REP1'; 
datosp2=REP2'; 

datos=[0 datospl (1) datosp1(2) 0datosp1(3) datosp1(4);0 datosp2(1) datosp2(2) 0 
datosp2(3) datosp2(4)] 

[Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos) 

%% Desplazamientos, giros 
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q= KT\Q 

%% Fuerzas en columnas 
Q2C=Q2(1:9,:); 

[FH]=fuerzas(ngl,ELEMC,LC,senC,cosC,VCC,E,q,Q2C) 

%% Fuerzas en vigas con cartelas 
Q2VC=Q2(10:21,:); 

[FC]=fuerzas_cartelas(ngl,ELEMVC,LVC,senVC,cosVC,VCVC,E,q,Q2VC) 



